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Abstract 

The effect of pressure on acceptor levels and hole scattering mechanisms 
in p-GaSe is investigated through Hall effect and resistivity measurements 
under quasi-hydrostatic conditions up to 4 GPa. The pressure dependence 
of the hole concentration is interpreted through a carrier statistics equation 
with a single (nitrogen) or double (tin) acceptor whose ionization energies 
decrease under pressure due to the dielectric constant increase. The pressure 
effect on the hole mobility is also accounted for by considering the pressure 
dependencies of both the phonon frequencies and the hole-phonon coupling 
constants involved in the scattering rates. 
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1 INTRODUCTION 

Gallium selenide (GaSe) belongs to III- VI layered semiconductor family which 
is characterized by a strong anisotropy in the chemical bonding. Its crystal structure 
consists of plane hexagonal lattices which are associated in pairs and may be 
stacked in different ways(/3, 7, e-polytypes)[jT[-0, being the most common one the 
corresponding to e-polytype {D\ h space group). The possibility of growing III- VI 
semiconductors thin films by van der Waals epitaxy |J- || opens new opportuni- 
ties to their potential practical applications, e.g. electronic devices^] and solar 
cells J7|. In this respect, the study of the fundamental electrical transport prop- 
erties is a very important issue. In the case of GaSe, at ambient pressure, these 
properties have been investigated in detail by different authors who have focused 



their attention on the role played by impurities such as CdpJ, Cu[H], N[1C], Sn[ll 



and Zn ||12|| . In contrast, very little is known about the behaviour of transport 
properties under pressure (P). 

Pressure experiments are an efficient tool to study III- VI semiconductors since 
they allow the tuning of the degree of anisotropy of these materials. Experiental 
studies of their optical(T^]-[|l^] and lattice dynamical JIT]] - pOf properties under 
pressure have been subjects of recent interest. However, systematic research of the 
electrical transport properties under compression has been only made for indium 
selenide (InSe), for which, resistivity (p) and Hall effect (HE) measurements under 
pressure have been reported very recently^!]]. 

In this article, we report a systematic study on transport properties under high 
pressure, up to 4 GPa at room temperature (RT), of the III- VI layered p-GaSe. In 
Section II we briefly describe the experimental setup. The data obtained on HE 
and resistivity measurements are presented in Section III. Finally, Section IV is 
devoted to the discussion of these results and their implications. 

2 EXPERIMENT 

The p-GaSe single crystals used in this study were grown by the conventional 
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Bridgman technique. Doping by tin (Sn) was performed by adding the pure element 
so as to get 0.05, 0.1, 0.2 and 0.5 at % in the stoichiometric melt of GaSe. Nitrogen 
(N) was introduced as GaN compound in a quantity to give a concentration of 0.5 
at % of N in the melt. We should point out that these concentrations are different 
from the effective ones, since during the Bridgman growth most of the initial 
impurity concentration is segregated to the end of the ingot J 



Samples with faces perpendicular to the c-axis were prepared from the ingots 
by cleaving and cutting with a razor blade. Thickness of the slabs was measured 
by using the interference fringe pattern in the infrared transmission spectrum. 
Typical sample dimensions were 70 /im in thickness and about 4x4 mm 2 in size. 
Pressure up to 4 GPa was applied to the samples by using the Bridgman cell that 



has been described in Reference [21]. In this case, we have used tungsten carbide 
anvils, 27 mm in diameter, which were put between the pistons of a 150-ton press. 
Sodium chloride was the pressure transmitting-medium. The pyrophillite gaskets 
were 0.5 mm in thickness, with a hole of 9 mm. The Bridgman gasket assembly 
and geometry are shown in the inset of Figure 1. Ohmic contacts were made in 
a van der Pauw geometry [f23|l by soldering silver wires with high-purity indium in 
gold contact pads which were previously vacuum evaporated. HE and resistivity 
measurements were carried out at RT. Direct electric current ranging from 50 /iA 
to 500 /iA was sent through the samples, and the resulting voltages were measured 
by a digital voltmeter. The linearity of the ohmic voltages on the injected current 
was checked out at different pressures. A magnetic field of 0.6 Tesla was applied 
parallel to the c-axis. The magnetic field was generated by means of a copper coil 
placed around one of the pistons of the press. 

3 RESULTS 

Figure 1 shows the pressure dependence of resistivity for N-doped and Sn-doped 
GaSe samples. The nominal doping concentration ([N] or [Sn]) of each measured 
sample is also given in this figure. Resistivity appears to decrease with pressure 
in all the samples here studied. This evolution is more pronounced in the samples 
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doped with N, in which p goes down a factor three in the explored range of pressure. 
This factor is nearly constant in all the Sn-doped samples, in which resistivity 
decreases with doping, except in the slabs with the highest tin content ([Sn] = 
0.5 %). 

The decrease of the resistivity with pressure turns out to be due to the increase 
of both the hole concentration and mobility, as shown in Figures 2 and 3. Figure 
2 shows the pressure behaviour of the hole concentration (jp) for different samples 
as determined through 

1 
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where q is the electron charge and Rh is the Hall coefficient. The Hall factor 



has been assumed to be Ifl0|, |1T] . It can be seen in Figure 2 that p non-linearly 
increases under compression, the relative variation of p being higher for samples 
doped with N. In samples doped with Sn, at zero pressure p is enlarged with 
increasing [Sn] with the exception of the 0.5% Sn-doped sample in which p is two 
orders of magnitude smaller. In addition, in this sample the relative change of p 
is the largest among the Sn-doped samples. 

As regards the mobility (/i), Figure 3 shows its pressure dependence. It can 
be seen there that the mobility also increases with pressure and that its value at 



ambient pressure is very similar to that obtained in previous works ||10|, [XT 

4 DISCUSSION 

4.1 Pressure Dependence of the Hole Concentration 

To understand the pressure dependence of p we have considered that N-doped 
and Sn-doped GaSe at RT are extrinsic as reflected in the temperature dependence 
of their transport properties ^ , |Tl| . Moreover, from the temperature behaviour 
of the hole concentration it was deduced that doping by N introduces only one 
acceptor level with ionization energy of E a = 210 meV|TD|. In the Sn-doped samples 



two impurity levels with ionization energies of E a = 155 and 310 meV have been 
observed [II]]. These levels appear to be connected with the presence of a double 



acceptor-impurity. This acceptor has been proposed to be an interstitial Sn atom 
in the octahedral interlayer site associated to a close Ga-vacancy, giving a local 
configuration similar to that of Sn in the layered compound SnSe2 ■ At RT the 
first level (E a = 155 meV) is not fully ionized, so that the pressure dependence of 
p is mainly determined by its behaviour under pressure. However, when increasing 
the doping of Sn to 0.5% a decrease of the hole concentration has been observed (see 
Figure 2) indicating that, at very high Sn concentration, a donor configuration of Sn 
in GaSe (isolated Sn interstitial or Sn substituting to Ga) becomes dominant and 
overcompensates the acceptor centers. Then, as in the heaviest doped samples the 
first level (155 meV) is completely compensated, the behaviour of the second level 
(310 meV) would determine the variation of p under compression. Therefore, we can 
analyze in both cases the pressure dependence of p through a single acceptor-single 
donor model for partially compensated p-type semiconductors [2^] . Within this 
model the hole concentration is given by: 
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where T is the absolute temperature, k is the Boltzmann constant, E a refers to 
the ionization energy of the acceptor level, Nd and N a are the donor and acceptor 
impurity concentration, respectively, and N v is the density of states of the valence 
band which can be written at RT as a function of the effective mass in the valence 
band m* = {ml^m^) 1 ^ 3 as: 

3/2 



N, = 2.509 x 10 



cm 



where m a , m* hLl and stand for the free electron mass, the perpendicular and 
parallel effective hole masses, respectively. Taking m* hi _ = 0.8 m and m^n = 0.2 
m at zero pressure^, pp] , we obtain m* = 0.5 m . 

In Eq.(2) one can see that p depends on the pressure through m*, E a and the 
term N d /2. Following the Kane model E7I], the variation of the effective mass with 



pressure is considered to be proportional to the change of the direct band-gap [|T5|]. 
For the pressure dependence of A^ it can be assumed that it is determined by 
the volume variation [0. Then, taking for N a , Nd and E a at ambient pressure the 



values deduced from the temperature dependence of the hole concentration |L0|, |TT 
and using Eq.(2), we can calculate E a at each pressure. The result is shown in 
Figure 4 (points). There it can be seen that the observed behaviour of p would 
be a consequence of a reduction of the ionization energies under pressure. For the 
same pressure variation (from to 4 GPa) the ionization energy of the N-related 
acceptor would decrease by 20 %, but that of the Sn-related one, only by 6.4 or 
8.1 %. 

The ionization energy of a single hydrogenic impurity level can be evaluated 
in an anisotropic crystal through the Gerlach-Pollman model By using the 



effective masses of GaSe from References [ER EHJ and the static dielectric constant 



at ambient conditions from Reference 17 one obtains a value of E a = 72 meV. 



The discrepancy between this value and the experiment can be attributed to the 
differences between the Coulomb potential used in the hydrogenic model and the 
true impurity potential. It has been shown that the true impurity potential can 
be represented by an effective potential which includes a central-cell correction to 
the Coulomb potential model fl2"9|. Then, to modelize the pressure dependence of 



E a we take the impurity potential to be[2£ 



1Z R* a* 

V(r) = 1 — exp(— br) + Br exp(— br) f(a,9) , (3) 

r L J 

where R*, a* and f(a,9) are the effective Rydberg energy, the effective Bohr radius 

and the anisotropy function defined in the Gerlach-Pollman model, respectively, 

b and B are the adjustable potential parameters, and Z is equal to 1(2) for the 

single (double) ionized acceptor-impurity. At small r, this potential looks like a 

sphericall well, with a depth Vo given by: 

n B + b 

V, = -q 2 ^— , (4) 

£o_i_£o|| 

where Eo± and £o|| are the static dielectric constant in the direction perpendicular 
and parallel to the c-axis, respectively. At large r it behaves like a screened 



Coulomb potential. The turning point is given approximately by 

B + b 



(5) 



Bb 

The proposed impurity potential model and the Coulomb potential are shown 
schematically in the inset of Figure 4. 

The Schrodinger equation of the system is solved by the variational method by 
using a trial function for the electron ground state (\l/o( r )) as: 

[¥ 

\E'o( r ) = V — exp(— for) . (6) 

V 71 

The variational calculation in which the total energy is minimized leads to: 

/ 5\ 2 16 
E a = -[-) R*Z%(a)--BR*a*Z (a) , (7) 



where Z (a) is the effective charge pEfl . Then the potential parameter B is chosen 



to fit the ionization energy of the single acceptor level at room pressure. The fitted 
value of B is shown in Table 1 together with Vq and r z . E a depends on the pressure 
through R* and a*. The first term of Eq.(7) decreases under pressure due to the 



reduction of R* as a consequence of the increase of £qii[117|, |30| . The second term 
also decreases under pressure, but more slowly because R*a* ~ (£ ||) -1 / 2 . 

The calculated pressure dependence of E a is shown in Figure 4 (solid lines). It 
can be seen there that the model fits quite well the experiment. The shift observed 
in the second Sn-related level may be connected with the fact that in the present 
model we have neglected the electron-electron interactions and, as a consequence 
of this, we could be overestimating the value of E a . To understand why E a suffers 
a greater decrease in N-doped samples, let us consider the parameters given in 
Table 1. It can be seen there that in this case the square well is deeper. This 
is just what one may expect because N is much more electronegative than Sn. 
That is why the central-cell correction is more important in N-doped samples than 
in Sn-doped samples. Because of this a deeper acceptor level is obtained. In 
addition, the second term of Eq.(7) would decrease more rapidly leading to the 
greater change observed in E a . 
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4.2 Pressure Dependence of the Hole Mobility 

Ionized impurity scattering and two-phonon scattering mechanisms must be 
considered in order to give a quantitative account of the pressure dependence of 
the hole mobility. The ionized impurity concentration has been assumed to be 



that obtained from the temperature dependence of the hole concentration |)T0|, |TT 
The phonons involved are the 138 cm -1 A' 2 homopolar optical mode and the 153.2 
cm -1 E' 3 LO polar optical mode[]I7|. The pressure dependence of the Frohlich 
constant 0] (a) and the hole-homopolar phonon coupling constant [|32| (g 2 ) were 
calculated from data found in the literature [0, |33], |3~4]| . As the hole scattering rate 
for LO phonons is obtained through an integration over all the possible directions 



of the phonon momentum, a was calculated through [|35 



2 1 

a = g«_L + ~a\\ . (8) 

We have obtained for both coupling constants a decreasing behaviour under pressure. 
The reduction of g 2 is basically a result of the increase of the homopolar phonon 



frequency under pressure [|17|, [Tgl - For a, in spite of the increase of an, due to 
the great increase of £o|| under compression [17, [3(|, the overall change is negative 
because of the decrease of a±, which is more than four times higher than an 



By means of an iterative method [24] we have calculated the pressure dependence 
of jji,. The ionized impurity scattering has been included through the Brooks- Herring 



relaxation time||36||, in which the ionized impurity concentrations has been assumed 
to be that used to calculate £„aa function of pressure. The results obtained for 
the samples doped with [Sn] = 0.05 % and 0.1 % and with [N] = 0.5 %, which are 
plotted in solid lines in Figure 3, agree quite well with the measurements. The 
isolated contribution of each scattering mechanism is also represented in Figure 3. 
The homopolar phonon scattering (curve 1) and the LO polar phonon scattering 
(curve 2) are the dominant mechanisms over the whole pressure range, but the 
ionized impurity scattering mechanism (curve 3) must be taken into account to 
reproduce quantitatively the absolute value of fi. The saturation of the increase of 
\x above 2.5 GPa is due to the saturation of the decrease of a as a consequence of 
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the compensation between the LO polar phonon frequency and Eq\\ increases. In 
addition, the increase of the concentration of ionized impurities, determined by the 
increase of p and the contraction of the volume under compression, leads to a re- 
duction of the impurity-limited mobility, which also collaborates to the saturation 
of \i. 

Finally, we want to point out that for samples doped with [Sn] = 0.2 % 
and 0.5 % the value of the zero pressure mobility cannot be reproduced with 
the present model. Samples from the 0.2 % ingot exhibit a phonon-controlled 
mobility in a larger temperature range than those from ingots with lower tin 
content. We think that this can be related to the fact that a large proportion 
of complex impurity centers are present in these samples, resulting in a reduction 
of the concentration of single ionized impurities JIT] . These complex centers ap- 



pear to have a lower scattering cross section than ionized impurities, producing 
that scattering by phonons is the dominant scattering mechanism in spite of the 
heavy doping concentration of the 0.2% Sn-doped sample. Nevertheless, the struc- 
ture of those centers and its influence on the hole mobility is not known and 
can hardly be included in our model. In samples from the 0.5 % Sn-doped ingot 
the compensation is very high and impurity scattering is dominant even at room 
temperature reducing the mobility to the low value observed (// = 9.5 cm 2 /Vs), 
but this value and its temperature dependence could not be accounted for through 
any known scattering mechanism. 

5 CONCLUSIONS 

Transport measurements have been carried out under pressure in N-doped and 
Sn-doped GaSe up to 4 GPa. Within the framework of a single acceptor-single 
donor model, the observed increase of p under compression has been interpreted 
to be due to the reduction of E a with pressure. Modeling the impurity potential 
we have obtained an expression for E a . This allows us to relate the decrease 
under pressure of E a to the increase of e ||- 111 addition, we have also discussed 
the different behaviour under pressure of the ionization energies of the acceptor 
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levels connected to Sn and N as doping impurities in GaSe. The higher reduction 
of E a observed in N-doped samples has been explained by means of the deeper 
central-cell correction obtained in this case. The observed increase of ji has been 
attributed to the decrease of the hole-phonon coupling constants. 
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Table 1: Model parameter B, deep of the square 
well Vq and turning point r z . 



Impurity 


B 


V 


r z 


Sni 


0.36 


-0.725 eV 


4.40 a* 


Sn 2 


0.36 


-1.120 eV 


2.47 a* 


N 


0.53 


-1.020 eV 


4.25 a* 
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Figure captions 

Figure 1: Resistivity as a function of pressure for different samples. The inicial 
doping concentrations in the growth solutions ([N] or [Sn]) are indicated in the 
figure. The inset shows the Bridgman gasket assembly. 

Figure 2: Hole concentration as a function of pressure for different samples: 
(• and |) [N] = 0.5 %, (V) [Sn] = 0.5 %, (A) [Sn] = 0.2 %, (□) [Sn] = 0.1 % and 
(o) [Sn] = 0.05 %. 

Figure 3: Pressure dependence of the hole mobility in several samples: (•) [N] 
= 0.5 %, (V) [Sn] = 0.5 %, (A) [Sn] = 0.2 %, (□) [Sn] = 0.1 % and (o) [Sn] = 
0.05 %. The solid lines are the results of our calculations. Curves 1 and 2 represent 
the homopolar phonon and LO polar phonon contributions, respectively. Curve 3 
represents the ionized impurity contribution as calculated for N-doped GaSe. 

Figure 4: Pressure dependence of E a . Curves labeled with N, Snx and Sn 2 
correspond to the nitrogen and tin levels, respectively. The solid lines illustrate 
the theoretically calculated dependence of E a according with the proposed model. 
The inset show a schematic comparison of the Coulomb potential and the proposed 
impurity potential model. 
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